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Abstract 

Different electron acceleration regimes in the evanescent field of a surface plasma wave are studied 
by considering the interaction of a test electron with the high-frequency electromagnetic field of a 
surface wave. The non-relativistic and relativistic limits are investigated. Simple scalings are found 
demonstrating the possibility to achieve an efficient conversion of the surface wave field energy into 
electron kinetic energy. This mechanism of electron acceleration can provide a high-frequency 
pulsed source of relativistic electrons with a well defined energy. In the relativistic limit, the most 
energetic electrons are obtained in the so-called electromagnetic regime for surface waves. In this 
regime the particles are accelerated to velocities larger than the wave phase velocity, mainly in the 
direction parallel to the plasma-vacuum interface. 

PACS numbers: 52.65.Cc,73.20.Mf,52.35.Mw 
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INTRODUCTION 


Electron acceleration by laser-plasma interaction has been studied extensively within the 
context of laser absorption by a plasma. It has also been studied for the development of 
techniques aimed at producing hot electrons in order to obtain improved energetic electron 
sources and to enhance secondary processes including X-ray, 7 -ray, positron production and 
ion acceleration [lj. In under-dense plasmas, various methods of wakefield acceleration [2H3j 
have been proposed and electrons up to ~ 4.2 GeV have been experimentally observed in the 
optimal configuration [5j. However, the electron acceleration mechanism invoked involves 
high laser intensity, well into the relativistic regime, and short pulse duration 77 = 40 fs 
resulting in a small current (total current of ~ 10 pC). It is thus of interest to search for 
alternative configurations wherein the current can be enhanced. 

In dense plasmas, an energetic electron population may be created by resonant ab¬ 
sorption |H]. In this case, the electrons are accelerated by the resonant plasma waves 
excited at densities around the critical density by the interaction of a “long” laser pulse 
with a gentle-gradient dense plasma. With the development of lasers of high intensity 
(JAq > 10 16 W cm _ 2 /im 2 ), short pulse duration (< lOOfs), and high contrast (~ 10 12 ), new 
electron heating mechanisms have also been proposed [MB- They involve very sharp den¬ 
sity gradients and over-dense plasmas which do not necessarily require a resonant plasma 
response. Among these mechanisms are vacuum heating [7j, and the so-called ponderomotive 
or J x B heating | 8 j. Numerical simulations [12j have shown that, in many cases, vacuum 
heating can be more efficient than resonant absorption. The advantage of these mechanisms 
is that they involve dense plasmas such that they bear the potential to generate very high 
currents. 

In order to improve the electron acceleration to relativistic values by a laser interacting 
with an over-dense plasma, particular targets designs can be used [13, [T4]. Also surface 
electron plasma waves resonantly excited by the laser on a structured target can be used 
to reach that goal [15, 16j. The electrons interact in this case with the surface plasma 
wave field that is highly localized at the vacuum-plasma interface and oscillates with the 
laser frequency. The local field amplitude is higher than that of the incoming laser such 
that the typical electron quiver motion is much faster than the electron thermal velocity. 
Experimental evidence for the feasibility to accelerate electrons and ions by this kind of 
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scheme have recently been reported mt 

Surface plasma waves can also be excited at the surface of a laser-produced ion channel 
(by a mechanism similar to wakefield acceleration). Their use has been proven to be an 
efficient means to improve fast electron generation 1. THIGH]- Moreover, full 2D simulations 
of J x B heating [21] have shown that these oscillations can be unstable and decay into 
a standing surface plasma wave, revealing at the same time the existence of a non-linear 
mechanism for generation of surface plasma waves. 

However, depending on the characteristics of the excited surface plasma waves the electron 
population can have very different features. In this paper we explore with a simple model 
the different acceleration regimes of an electron in the evanescent electromagnetic field of a 
surface plasma wave by considering the interaction of a test electron with the high-frequency 
surface wave fields. The non-relativistic and relativistic limits are investigated by means of 
ID and 2D test particle simulations in order to find the optimal regime for an efficient 
conversion of the surface wave field energy into electron kinetic energy. The ID dynamics 
calculations are mainly performed in order to analyze the motion of the electron in the 
direction perpendicular to the surface and its excursion into the evanescent surface plasma 
field. Some analysis of the ID motion in the direction parallel to the plasma-vacuum surface 
is also performed. 

After a brief summary of the structure of the surface plasma waves (section II), the 
electron acceleration mechanism proposed is studied, first of all for non-relativistic over- 
dense plasmas (section III). The treatment of this case is based on classical concepts such 
as the idea of a ponderomotive potential [22] and the separation of high-frequency and low- 
frequency effects, which allow an analytical ID treatment (subsection III-A). This part is 
then complemented by full 2D test particle simulations in subsection III-B. It permits to 
identify two different situations for the electron acceleration: the electromagnetic regime and 
the electrostatic regime. This non relativistic part of the study is also a tool to understand 
the results in the relativistic regime. The next section (section IV) studies the relativistic 
limit where the quiver velocity becomes close to the velocity of light. First a numerical 
ID study is presented (subsection IV-A), as within this domain the possibilities for an 
analytical treatment are limited. This allows the determination of the different regimes, 
and of the optimal conditions in terms of electron acceleration. To complete the study, a 
full 2D numerical study is presented in subsection IV-B. We find that the electromagnetic 
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regime is the most efficient and taht an electron initially at rest can be self-injected and 
phase-locks on the vacuum plasma side through vx B mechanisms. This way the electron is 
accelerated to velocities larger than the wave phase velocity, mainly in the direction parallel 
to the plasma-vacuum interface. The resulting electron total energy scales as 7 v a sw mc 2 with 
a sw = eE sw /mcu where E sw is the SPW field component perpendicular to the surface of the 
plasma, = 1 / \/l — (u f /c) 2 with v v the phase velocity of the SPW and u is the surface 
plasma wave frequency. 


STRUCTURE OF THE SURFACE PLASMA WAVES 

We consider hereafter a 2D homogeneous plasma in the (x, y ) plane, which supports a 
surface plasma wave that propagates at the plasma-vacuum interface (along the y-direction). 
The vacuum is along x < 0. On the vacuum side (x < 0), the resonant surface wave field 
has the form: 

/ 

E x>v = —E sw f (t, x)sin(ut - k y y + 0), 

< E VjV = E sw jr-^f(t, x)cos(u)t - k y y + 0), (1) 

-^ 2 ,v E sw -^f(t, x)sin(ujt k y y T 0), 

with: 

f(t,x) = exp(x/ Le )V ) exp(— 2 t 2 /r^ w ), 

while on the plasma side (x > 0 ): 

/ 

E sw x)sin{u)t k v y T </>), 

^ E sw ~ 9(^1 x)cos(ujt k y y T 0), (2) 

B Z , P = E sw v -fg(t, x)sin(ujt - k y y + 0), 

with: 

g(t,x) = exp(-x/L EtP ) exp(- 2 t 2 /t 2 J, 


where E sw is the maximum value of the electric field on the vacuum side in the ^-direction. 
In the following this field will serve as reference field. t sw is the mode lifetime, 0 is the 
phase, v p = uj/k y is the phase velocity of the wave and Lp j V and Lp, p are the evanescence 
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length in the vacuum and plasma respectively. The expression for the evanescence length 


in vacuum is given by L E , V = y (k 2 — cn 2 /c 2 ) -1 . The evanescence length within the plasma 
is typically significantly smaller than the evanescence length in vacuum [23j. It is given by 


Le, p = y (ky — (co 2 — cu 2 pe )/c 2 )-' < L EtV . t sw will depend on the damping mechanisms that 
affects the plasma wave, such as linear and non-linear Landau damping, and wave breaking 
effects. 

In the cold-plasma limit, the thermal corrections to the dispersion relation for the surface 
plasma waves can be neglected, and in the non-relativistic limit we have [23J the relation: 

1 (id/uj pe ) 


k 2 y C? 


( 3 ) 


UJ 2 1 — 2{uj/uo pe ) 2 ’ 

where u> pe = Anne 2 jm is the electron plasma frequency, n the electron density and m the 
electron mass. By solving Eq. ^ for cu(k) it is easily seen that there is an upper limit for 
uj < oj pe /\/2. The surface plasma waves which satisfy this dispersion relation are localized 
at the plasma surface. 

Let us now consider in the following some limiting situations for the surface plasma 
wave. The relevant parameter to identify the different regimes for the electron plasma wave 
is u)/cu pe . We will therefore consider different values of this parameter, and the corresponding 
wave vector that satisfies the dispersion relation Eq.Q. 

In the limit uj <C uj pe , which is the so-called electromagnetic limit for the surface plasma 
waves, the magnetic field is of the same order of magnitude as the perpendicular electric 
field ( E x ), while the component of the electric field parallel to the surface (E y ) is negligible 


(see EqjlJ which in this limit yields \E y 


^f-\E x \). Thus the electric field is mainly in the 
x direction (i.e. perpendicular to the plasma-vacuum interface). It should be noted that in 
this case two different situations can be expected, depending on the field intensity. At low 
field intensity, it can be anticipated that the role of the magnetic field will be negligible, 
such that a ID model can be used, wherein the motion only takes place along x. In the 
opposite case that the value of the fields is large enough to accelerate the electrons into 
the relativistic regime, the magnetic field contribution will play an important role. In the 
electromagnetic limit where uj <C uo pe , the wave phase velocity along the surface, v v = uj/k !n 
is less than but of the same order of magnitude as the velocity of light. Moreover the 
evanescence length of the wave in vacuum can be quite large: the smaller the value uo/cu pe , 
the larger the vacuum evanescence length. For example, if uj/uo pe = 0.05 we have from Eq. 


5 






TABLE I: Summary of the parameters for the surface plasma wave with the following definitions: 


K = ( kyc/oo ) 2 , Le, v = w 2 /c 2 ) 1 , Vy = c/\fK and 7 ^ = \/The different 

values of uj/cu pe correspond to different regimes for the surface waves. 


u/uj pe 

K 

Le,v/^ 0 

-^y,v / -^x,v 

v v /c 

W 

0.7 

25.5 

0.032 

0.98 

0.198 

1.02 

0.6 

2.28 

0.14 

0.75 

0.662 

1.33 

0.22 

1.053 

0.69 

0.2 

0.974 

4.54 

0.05 

1.0025 

3.18 

0.05 

0.998 

20 


|3j) K = kyC 2 /u > 2 = 1.0025, and Le, v = A 0 /27 \\[K — 1 = 3.18A 0 , for A 0 = 27 tc/uj. If instead 
we consider uj/oj pe = 0.22, we have K = 1.053, and L EfV — 0.7Ao- As we will see in the 
following, large phase velocities v.„ = c/s/K , and large evanescent lengths L E}V are favorable 
situations for efficient electron acceleration within the highly relativistic regime. Note that 
in the electromagnetic limit the evanescence length in the plasma is much less than the 
evanescence length in vacuum. In particular in the limit u> <C u> pe we have L E . p ~ c/co pe . 

On the other hand, in the limit uj ~ u> pe /^/2, we are close to the electrostatic limit where 
the magnetic held is always negligible. The two components of the electric held are of the 
same order of magnitude and the evanescent length is very small. Typically we always 
have L Ev ,L Ep <C Ao- In particular, if cu = 0.7 u> pe ~ u> pe /s/2 , we have K = 25.5, and 
L e , v = 0.032A 0 . In this case, the wave phase velocity along the surface is much smaller than 
the velocity of light. 

We have thus identihed two important parameters characterizing the surface plasma 
waves, the phase velocity, v v , along the surface and the evanescent length L E v perpendicular 
to the surface. These two parameters will allow us to dehne different regimes governing the 
particle acceleration. This will be done in the next section, where we will concentrate on the 
solution of the wave on the vacuum side, since it corresponds to a more favorable situation 
for electron acceleration. In table I, we have reported these parameters for different values 
of u/co pe to be discussed below. 


6 












NON-RELATIVISTIC LIMIT 


We first analyse the case of low surface plasma field intensity such that a sw = 
eE sw /mcjjj <C 1. Relativistic effects are then negligible. Let us recall here that E sw is 
the maximum value of the electric SPW field component perpendicular to the surface in 
vacuum, E XtV , and that the perpendicular component of the field is significantly reduced in¬ 
side the plasma. For this reason we focus on the motion of the electrons in the vacuum side. 
As we are typically in the limit t sw 2i t/lu we will have a slow time evolution combined 
with a fast time variation due to the high-frequency collective electron oscillations. 

Acceleration perpendicular to the surface 

We start the study of the motion of the electrons by considering only the behavior in the 
perpendicular direction, x, with the aim of highlighting the role of the evanescence length 
of the field. In this case, we consider a fixed value of y , viz. y = 0, and we neglect the 
contributions of E VtV and B zv . This is quite good an approximation for the electromagnetic 
case since E VjV /E XtV <C 1 (see table I) and the electron velocity is much smaller than the 
velocity of fight. Associated to the maximum field amplitude we can define the electron 
quiver velocity v osc = eE sw /'y osc mu>, where 7 0SC = -y/l + (eA^/cmcu) 2 . In the non-relativistic 
limit considered here 7 0SC ~ 1, but we introduce it here in order to formulate the most general 
definition of v osc . 

The high-frequency motion of a single electron in the electric field will be characterized 
by the length A = v osc /lo. This quantity should be compared with the length scale at which 
the field L E , V varies and we define hereafter the ratio Rl = L E , V / A. 

If Rl 1, that is when v osc is not too large, the electrons will have time to perform 
many oscillations before leaving the resonant field. This is the relevant physical regime at 
low laser intensities where relativistic effects are negligible. In this case the exact values of 
L E:V , or uj / (jj pe do not matter, since to a good approximation all frequencies will correspond 
to the same physical situation, as long as Rl 1. In this limit, the electrons undergo 
the effect of the ponderomotive force, which implies an averaging over the fast oscillatory 
motion in the field gradient. Hence, the typical kinetic energy acquired by an electron can 
be equal to the ponderomotive potential U osc = (1/4 )m e v‘l sc . 
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In fact, the lifetime r sw of the surface plasma wave may play an important role in de¬ 
termining the amount of net kinetic energy that can be gained by the electrons within the 
surface plasma wave. It has been shown in particular [ 23 ], that only a partial conversion of 
the potential energy into kinetic energy can occur if t sw is shorter than the time (~ L EtV /v osc ) 
needed by the electron to explore the whole spatial extension of the surface plasma wave 
held. In this case, the net kinetic energy gained by the electron will be only a fraction of 

UqSC ■ 

There exists an additional source for electron acceleration, that is not of ponderomotive 
origin inasmuch it does not specifically involve a low-frequency time scale or a spatial depen¬ 
dence of the held. At the moment they enter into the electric held of the surface wave[24], 
the electrons can also gain some kinetic energy on a sub-period time scale, that is, before 
even feeling the effect of the ponderomotive force. This extra kinetic energy may be simply 
related to the motion of the electron entering in a spatially constant oscillatory held and 
is, therefore, a strong function of the entry phase of the electron in the surface wave. The 
variations in the entry phase rehect the variations of the entry times of the electron into the 
surface wave held. However in the following we keep the initial time constant equal to zero 
and we consider explicitly different values of the phase 0. 

A hrst value for the total kinetic energy acquired by an electron can be obtained from the 
zero order term in a series expansion for the solution of the equations in the presence of an 
electric held gradient. Thus we can solve the equations of motion of an electron with initial 
velocity v = —v 0 at t = 0 in a spatially constant held, whereby v 0 is of the order of the 
thermal velocity. After averaging over the fast motion, the resulting velocity as a function 
of the phase is given by: 

< v{4>) >h.f= -v 0 + v osc cos 4>. (4) 

Formally this is equivalent to having as initial velocity < v(ip) >h.f. instead of — v 0 . The hrst 
order solution in the electric-held gradient, yields the ponderomotive force: At this stage 
any other dependence on the phase disappears in the averaging procedure that is necessary 
to dehne the ponderomotive energy, such that the hnal kinetic energy (Wkj) of the electron 
is obtained by adding the ponderomotive potential to the kinetic energy associated with the 
velocity given above. By assuming that the ponderomotive potential energy is completely 


converted into kinetic energy, we obtain: 


Wkj — U osc + 2 U a 


v 0 


COS 0 


( 5 ) 


It is convenient to express Eq. ([5]) in terms of the final value Vf of the velocity acquired by 
the electrons: 

/ N 2 \ V 2 


°- 5 + r^ +cos0 ) • (6) 

This equation shows that the phase for which the electron obtains the maximum energy 
(best phase), is given by 0 = ir. For this phase, the electron leaves the surface plasma wave 
field with a maximum kinetic energy such that Vf is greater than v osc /\/2. Hence, due to 
this phase dependence, there will exist electrons which are accelerated to energies higher 
than the ponderomotive potential energy. Even if the contribution to the electron energy 
in this case is not uniquely of ponderomotive origin as explained above, we call this the 
ponderomotive regime. 

In the following, we have numerically solved the ID equation of motion: 


+ < 7 > 

with = jrnRI for an electron subjected to an external field ~S, representing the field of a 
surface plasma wave, as defined in section I. The initial conditions are such that at t = 0, 
x = 0,u = — vo, for vq corresponding to lOOeF, that corresponds to vo/v osc = 0.16. The 
surface plasma wave lifetime is taken as 30r 0 where r 0 = 2ir/u>. The electron motion is solved 
by using the Vay pusher [26]. We have adopted a value E sw = 2.76 x 10 9 V/cm (a sw ~ 0.086), 
for which 7 osc ~ 1. The numerical results are reported in fig.l for u)/cu pe = 0.22. For this 
case we have R E ~ 50 and L E)V = 0.7pm. For comparison, the analytical curve obtained 
from Eq. ([6j) is also reproduced in the figure. 

For the phase 0 ~ n, the numerical and theoretical results are in very good agreement. 
They correspond to electrons that have acquired the possibility to fully explore the field gra¬ 
dient during their motion, after having been accelerated to their maximum energy. When 
the initial phase approaches 37 t/ 2 , the discrepancy of the numerical curve with the theoret¬ 
ical prediction of Eq. 0 grows larger. This prediction is obtained by assuming that the 
ponderomotive potential is completely converted into kinetic energy. This can be explained 
by noticing that for the less favorable phase the electrons are moving slower such that they 
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FIG. 1 : Non-relativistic limit: final electron velocity Vf/v osc as a function of the entry phase (f) in 
the SPW field. Here \mvQ = lOOeV and v osc /c ~ a sw = 0.086. We plot equation (6) (black full 
line), and the numerical values for Rl ~ 50 = 0.22 u pe ) (*), and Rj j ~ 10 (u = 0.6cj pe )(*). Only 

the electrons in the vacuum side have been considered. 

need more time to cross the field. Under these conditions, the fact that the wave has a 
finite lifetime can no longer be neglected ( ujt sw ~ 180) such that only a partial conversion 
of potential energy into kinetic energy can take place. This is verified by considering an in¬ 
termediate case, with a shorter evanescence length Le, v '- for example Rl ~ 10 ( u> = 0.6 u> pe , 
Le,v = 0.14 fim). As expected, the acceleration obtained and reported in fig|l] is now closer 
to the theoretical prediction for ~ 37 t/ 2. In conclusion, we can say that the phase of the 
field experienced by the electron plays an essential part in determining the range of energy 
that can be acquired on traveling through the field. Moreover the role of the surface plasma 
wave lifetime is determined by the parameter L E:V /v osc t sw , and is negligible if this parameter 
is small. 

2D case 

We now consider the full 2D motion of the electron in the non-relativistic limit in more 
detail. Note that in this case, the entry phase 4> is both representative of electrons entering 
the wave at different times or in different space locations in the ^-direction. (In the simu¬ 
lation, only the electrons in the vacuum side have considered and we have taken x = 0 and 
y — 0 at t — 0). 

In the electromagnetic regime where cu <C co pe , taking into account the E y and B z compo- 


10 


































FIG. 2: Comparison of the final electron velocity Vf/v osc from the ID and 2D simulations as a 
function of the phase <p at entry of the surface plasma wave field for Rl ~ 50 (u = 0.22u pe ) (upper 
part) and Rl ~ 10 (u = 0.6cj pe ) (lower part), o ID simulation, ♦ x component of the momentum 
and * total momentum in the 2D simulation. Only the electrons in the vacuum side have been 
considered. 

nents of the field is expected to have a minor effect on the particle motion for two reasons: i) 
The particle speed is always much smaller than the speed of light, c, and thus the magnetic- 
field contribution is small, and ii) E V:V /E XtV <C 1. This can be observed in the upper part of 
fig|2] where for comparison we have reported the final electron velocity from the ID and 2D 
simulations as a function of the phase <fi when the electron enters the surface plasma wave 
field for Rl ~ 50 (cu = 0.22 u pe ). In this case E y:u /E xy] = 0.2. We can see that there is a 
very small difference between the ID limit of the momentum and the x component of the 
momentum in the 2D limit. As in the one-dimensional case, we have a well-defined shape for 
the velocity distribution as a function of the phase. The existence of a flat extremum around 
<f> ~ 7T implies that many particles will have the same value of the momentum despite their 
difference in phase. This results in a bunching of the electrons in the direction perpendicular 
to the plasma surface in momentum space, and a concomitant energy bunching [2.5J. 

Increasing co/Lu pe we enter the electrostatic regime where 2D effects may become impor¬ 
tant. For Rl ~ 10 (uj = 0.6co pe ~ uj pe /\/2 ) we have E y ., } / E XjV = 0.75. A significant difference 
between the velocity in the ID case and the x component of the velocity in the 2D case can 
now be observed as displayed in the lower part of fig j2j The analysis of the ratio p y /p x of the 
final electron momentum for the 2D case as a function of the entry phase (f) (see fig j3] (right)) 
shows that the weight of p y is especially important for 4 < <fi < 5. The ID description is 
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FIG. 3: Evolution of the ratio p y /p x of the final electron momentum component for the 2D case as 
a function of the phase </> at entry of the surface plasma wave field for Rl ~ 50 {uo = Q.22uj pe ) (left) 
and Rl ~ 10 (cj = 0.6uj pe ) (right). Only the electrons in the vacuum side have been considered. 

then no longer valid since p y becomes of the same order of magnitude as p x . With respect 
to the estimate based on the ponderomotive potential Eq. a, the overall velocity acquired 
by the electron is enhanced by an amount that varies between 15% and 60% (depending on 
the phase). In the left part of fig|3j we have for comparison also plotted the ratio p y /p x in 
the electromagnetic limit (a; = 0.22cv pe , Rl ~ 50). It can be observed that in this limit p y 
is always significantly smaller than p x . We point out that in the case Rl ~ 10, the final 
momentum acquired by the electron is higher than in the ID case, even if the final p y is 
close to zero for phases between 3 and 3.5. This is a result of the history of the particle 
motion as at earlier times p y exhibits an oscillatory behavior that influences also the motion 
in the x direction. 

To summarize the case of low surface plasma field intensity, the electron motion within 
the electromagnetic regime is essentially ID in the direction perpendicular to the surface. 
If the SPW lifetime is higher than Z/£; ; „/u osc , the electrons with a favorable initial phase will 
gain more energy than can be accounted for on the basis of the ponderomotive contribution 
only. In the electrostatic regime, the energy gain will be even higher (~ 25% increase) due 
to 2D effects. 

RELATIVISTIC REGIME 

We consider now very high field intensities such that a sw = eE sw /m<xo > 1. The electrons 
can then be accelerated to relativistic velocities. (Here E sw is the maximum value of the 
electric field in the x-direction in vacuum and we have taken it for reference field). The 
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electrons will then have a behavior in the field that varies, depending on the issue if for 
a sw 1 the field evanescence length Le. v becomes comparable or shorter than the charac¬ 
teristic high-frequency electron motion length A = v osc /u> ~ c/uj. In this limit various other 
parameters become important such as the wave phase velocity, v v , and the relative value of 
the two components of the electric wave field at the surface, E yv /E xyu . 

In this section we first present a parametric study (section IV. A) of the energy transfer to 
the electrons. Only the transverse field of the evanescent surface wave is thereby considered. 
Within this ID approach it is possible to make some analytical estimates which allows 
identifying the different interaction regimes and the role of the evanescence length. The 
full 2D simulations follow in section IV.B, where the importance of the acceleration in the 
direction of the propagation of the surface wave will become evident. 


Acceleration perpendicular to the surface: role of the evanescent length 


Due to the high values of the surface wave field and the very short life time of the 
modes under consideration, we have here v osc ty, where iq might be in the range of 
~ lOOeV — 1 keV. We therefore do not expect our results to have any significant dependence 
on the initial value of the electron velocity vq, which is therefore neglected in the following 
analytical estimate. 

In the relativistic electrostatic limit where uo ~ c o pe /\[2, we can give an analytical estimate 
of the momentum acquired by the electrons as a function of the entry phase, since Rl <C 1. 
In this case, an electron can acquire a large velocity v < v osc ~ c over a time span less 
than a period, such that A Le. v - This means that the electron will leave the field 
and move towards the vacuum before the field will have had the time of going through an 
oscillation. The maximum energy gained by the electron under this hypothesis is given by 
Wkj = eE{4>)LE,vi since we can assume that the electron sees a roughly constant value of 
the electric field during the time it spends in the field of the surface plasma wave. 

Consequently, the value of the normalized electron momentum Pf/mv osc in this regime is 
given by (for 7r < < 2n): 


Vj_ 

me 



Pose Le^v 

me Aq 


cos {<j> + |) 



( 8 ) 
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The validity of this formula has been checked numerically for the case Rl = 0.2 (oj/cu pe = 
0.7), for which the hypotheses leading to the equation above apply. In figjd] we plot pf/mc 
versus the entry phase for Rl = 0.2 ( uj = 0.7uj pe , Le = 0.032 pm, a sw = eE sw /m<xo = 8.6) 
and a surface plasma wave lifetime t sw of 3 to, where To = 2i t/oj is the period of the wave. 
Notice that the result presented in the figure is in fact basically independent of the SPW 
lifetime, since the electron spends less than a period in the SPW held. As we see, the 
analytical formula Eq.(|8]) can be considered as a fairly good estimate of pf versus </>. In 
particular it provides with very good accuracy the maximum value of pf that can be obtained 
by an electron in this regime characterized by Rl <C 1, and thus the lower limit for particle 
acceleration at a given wave amplitude compared to optimum 2D situation. In that sense it 
can be used as a reference case, also owing to the fact that the energy transfer has a simple 
analytical formula. 



FIG. 4: Final electron momentum pf/(me) as a function of the entry phase in the surface plasma 
wave field for Rl = 0.2 (u = 0.7w pe )(the black ful line is the analytic prediction Eq|8]). Only the 
electrons in the vacuum side have been considered. 

When the scale of variation of the held is of the same order or larger than the distance 
explored by the electron in the held during one period, we have Rl 1. We denote 
in the following this situation as the relativistic electromagnetic regime. In this case, the 
kinetic energy acquired by an electron depends strongly on the phase of the held. As the 
possibilities of an analytical treatment are limited here, we present a numerical investigation 
of the hnal value of the momentum as a function of the entry phase. In the following we 
consider values of Rl ranging from 0.9 to 20, and we solve the relativistic equation of motion 
of electrons entering in the oscillatory evanescent held of the SPW, as dehned in section III. 
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First of all, the assumption made that in this relativistic regime the details of the initial 
distribution function of the electrons do not affect the final range of momentum acquired by 
the electrons leaving the field, has been checked. This assumption is already well verified 
for a sw = 2.72 and T e = 100 eU — 1 keV. Consequently, an uncertainty about the initial 
temperature of the electrons at the plasma surface will not significantly affect their final 
energy. In the following we take thus a constant value of vq corresponding to lOOeU. 

In this section we have taken electrons entering the field during one period only, at 
the moment that the time envelope of the surface wave is at its peak (t=0). Under this 
hypothesis the only effect of the surface plasma wave lifetime is reducing the amplitude of 
the field seen by the electrons during their motion towards the vacuum region, as was already 
observed in the non-relativistic case. As v osc is now close to the velocity of light, we have 
in most cases L e ^ v /v osc t sw -C 1, such that we can neglect the effect of the finite lifetime of 
the surface plasma wave in the discussion. In all the simulations reproduced here its value 
is taken equal to r sw = 30to. 

We consider hereafter first the strong surface wave field limit with a sw = 8.6 and values 
of Rl ranging from 0.9 to 4.3 (fig{5]). In a second stage we consider values from 4.3 to 
20 (figj 6 |) in order to fully explore the relativistic electromagnetic regime. Notice that in 
the relativistic limit the quantity Rl = L E v co/v osc , depends only on the ratio co/cu pe . In 
fact, since v osc — > c, we have, using the surface plasma wave dispersion relation, eq. (|3]) 
R l ~ L e , v u/c = co pe /co^Jl - 2 cn 2 /u; 2 e = 7 ^/VK ~ 7 

As it can be seen from figj 5 j the maximum value of the momentum pf , that can be 
acquired by an electron, increases steadily with R L . The maximum value is obtained for 
Rl = 4.3, for which we have L E v < A ~ c/co, and Pf max > p osc ■ It can also be seen that an 
electron with an entry phase close to or larger than 37t/2 will acquire a larger momentum 
than those with neighboring values for the phase, resulting in a minimum, and then a local 
maximum. This can be related to the fact that such electrons make a turn after traveling in 
the field gradient, such that instead of being pushed back into the plasma they start moving 
away towards the-low field region. When R L is close to one, this acts like a sort of resetting 
of the initial conditions, as illustrated in ref. [25]. 

Let us now consider the results from fig| 6 j where we have plotted pf/(mc ) versus the 
entry phase for larger values of Rl- It is clear from the figure that the maximum value of 
Pf does not change much as Rl increases. The maximum energy transferred to an electron 
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FIG. 5: Final electron momentum pf/(mc) as a function of the entry phase in the surface plasma 
wave field. Various values of uj/uj pe are considered. By varying this parameter we explore the 
values from Rl — 4.3 (uj/uo pe = 0.22) to Rl = 0.9 (uo/uj pe = 0.6), where the latter value tends 
towards the electrostatic regime. Only the electrons in the vacuum side have been considered. 
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FIG. 6: Final electron momentum pf/(me) as a function of the entry phase in the surface plasma 
wave field going from Rl = 4.3 (uo/uo pe = 0.22) to Rl = 20 (w/uj pe — 0.05). Only the electrons in 
the vacuum side have been considered. 

diminishes slightly with respect to the case where Rl > 1, ^ c/cc, for which pf > p osc , 

but as Rl gets larger we still have pf ~ p osc . The “hills” that appear in the figure for 0 
close to 37t/2 for the case Rl — 10 are again due to particles that perform a turn in vacuum 
before leaving the high-field region, a single turn for particles in the first “hill” or two turns 
for the particles that form the second “hill”. However we see that there are no such “hills” 
if Rl — 20. This is because now the scale of variation of the field is much larger than the 
path length covered by an electron over a single period, such that the particles either are 
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pushed back into the plasma or perform many turns before leaving the high-field region. 
In this case the effect of “resetting” the initial conditions which is able to produce a local 
increment of pf around <j> = 3tt/2, does no longer show up. 

To summarize: In the relativistic regime, the limit Rl > 1, L Ev ~ c/u corresponds to a 
much more efficient mechanism to convert the energy from the plasma wave to the electrons 
than in the opposite limit Rl ■ Cl. We have not reported here the results for finer variations 
of u)/ijOp e , but as long as Le, v is close to ~ c/u, the maximum value of pf does not change 
significantly. Moreover, the condition Rl > 1 is necessary to obtain an efficient acceleration 
of the particle along the surface, as this condition permits the particle to stay for a long 
time close to the surface and to interact with the field parallel to the surface. This will be 
shown in the next section. 

2D simulations 

The quantity E y ^ v /E XtV , reflects the relative weights of the field components in the parallel 
y and the perpendicular x directions. From an inspection of its values reported in table I, we 
expect that the influence of the parallel component of the field on the electron motion will 
vary according to the regime considered. In this section, we are thus going to investigate 
the full 2D motion of the particle. The simulations have been run taking the finite lifetime 
of the surface plasma wave equal to t sw = 50to ~ 314a; -1 . 

We first consider the relativistic electrostatic case, corresponding to Rl -C 1 and 7 ^ ~ 1. 
As found in the previous section, this is the case for which the energy transfer mechanism 
is less efficient as the particle spends only a very short time in the perpendicular field. 
However, the two-dimensional effects lead to a final total momentum that is larger than 
in the case where only the transverse field is considered. This can be observed in fig{7^ 
where we reproduce the total final electron momentum for the case Rl = 0.9 (u/u pe = 0.6 
and 0 V = 1.33). This figure also displays the value of the component of the final electron 
momentum along the perpendicular direction and the value for the momentum obtained 
from a ID model. Like in the non-relativistic regime, the effect of the electric field parallel 
to the surface is not negligible, especially for values of <f> ~ 4. This can be seen from table I. 
The effect is of the same order magnitude as that of the field in the perpendicular direction. 

The maximum value of the momentum in the 2 D case is roughly twice as large as the 
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value obtained in the previous section. An inspection of the ratio p y /p x reveals that the 
motion parallel to the surface is no longer negligible in the case where the electron energy 
gain is maximal. Although this motion contributes to the enhancement of the maximal 
energy reached, the energy gain in this regime remains of the same order of magnitude as 
in ID. 
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FIG. 7: Comparison of the final electron momentum for a sw = 8.6 from the ID and 2D simulations 
as a function of the entry phase </> in the surface plasma wave field a) for Rl = 0.9 ( 7 ^ = 1.33, 
Lo/uipe = 0.6) b) for Rl = 4.3 ( 7 ^ = 4.54, Lo/u pe = 0.22) and c)for Rl = 20 ( 7 ^ = 20, uj/uj pe = 0.05). 
o ID simulation, ♦ x compoment of the momentum and * total momentum in the 2D simulation. 
Only the electrons in the vacuum side have been considered. 


We next consider the relativistic electromagnetic regime wherein R L and 7 ^ increase from 
values > 1 to values W 1 . As in the previous section (IV. A) we first take Rl = 4.3 
( 7 p = 4.54 and co/oj pe = 0 . 22 ). As can be observed in fig{7j) the final value for the maximal 
electron momentum is significantly larger than the ID value. From the fact that the p x 
component in the 2D case is very close to value of the momentum in the ID case (p x ~ p osc ) 
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we can infer that most of the acceleration for the optimal entry phase 0 = 3.4 takes place in 
the parallel direction. This is corroborated by fig{8^ : the trajectory of an electron for the 
phase 0 = 3.4 from time cut = 0. to cut = 40 is reproduced superposed to the E y field of the 
surface plasma wave at the final time. In this case the electron always sees an accelerating 
field (multimedia view). The trajectory of an electron in the same time interval for the 
lowest energy entry phase 0 = 5. is shown in fig j8]b. In this case the electron sees both 
an accelerating and decelerating field, acquiring mainly perpendicular energy (multimedia 
view). A similar behavior is observed for the case Rl = 20 in fig{7}-. However, here the 
maximum energy is even larger (mainly due to an increase of the parallel momentum as we 
still have p x ~ p osc ). 


0 

a) 


15 
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30 

-15 



b) 



FIG. 8: Trajectory of an electron in the surface plasma wave field from time uot — 0. to uot — 40 
(before it leaves the accelerating region). In a) the optimal entry phase (j) = 3.4, that always sees 
an accelerating field (Multimedia view ); in b) the lowest energy entry phase </> = 5.0, that sees 
both an accelerating and decelerating field, acquiring mainly perpendicular energy (Multimedia 
view ). Parameters correspond to fig. [tJd. The electric field E y at time uot — 40 is superimposed 
to the trajectories for reference. 

In order to accelerate an electron in the direction parallel to the surface by a wave that 
propagates along the surface it is necessary that the particle velocity be close to the wave 
phase velocity. This can be obtained by injecting the electron right at the beginning with 
the appropriate velocity. Alternatively an electron may be initially at rest, and then become 
self-injected and phase-locked thanks to some other mechanisms. In our case, a significant 
contribution to the self-injection comes from the v x x B force. This is shown in figj9] where 
we reproduce the v x and v y components of the electron velocity as a function of time for 
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the phase, cp ~ 3.1, which gives the higher p/m e c value in fig {?]::. For comparison we also 
plot the v y component for a case where E x and B~ are set to zero, such that the force is 
absent. While in the latter case the electron oscillates in the field without gaining energy 
on average, we can see that the particle gets an extra contribution in the y direction and 
that its velocity becomes close to v v ~ c when the full 2D calculation is included. 



FIG. 9: Evolution of the normalized electron velocity with time for the entry phase ^ = 3.14 in 
the surface plasma wave field within the relativistic electromagnetic regime for Rl = 20 (yy = 20 
and Lu/u)p e = 0.05): v x , v y (2D simulation) and v y (setting E x = B z = 0). Only the electrons in 
the vacuum side have been considered. 

To further analyze the electron behavior in the SPW field, it is now useful to study the 
particle-field interaction in the reference frame of the wave phase velocity. To this extent, 
we perform a Lorentz transformation to the frame of the surface wave moving at iy = uj/ k y . 
In this frame, we have: 

/ 

cut - k y y = -k y y / 7 V , 
i= 7 <p(t - v v y/ c 2 ), 

<7 = llv( 1 -v< P Vy/c 2 ), (9) 

X — X, 

z — z, 

where the tilde refers to the quantities in the moving frame. When we neglect the lifetime 
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of the surface plasma wave, the field expressions on the vacuum side Eqs. (JT[) become: 


E x , v = -E sw ^-sin(-k y y/-/ lf + 0), 

< Ey )V = E sw cos(—k y y/ 7 ^ + 0), (10) 

B z , v = 0 , 

\ 

with h(x) = exp (x/L E}V ) with x < 0 . 


It should be noted that after this Lorentz transformation we now are dealing with 
an electrostatic problem on the vacuum side in the moving frame, with E XjV ~ E lhV . The 
electrostatic potential, U, defined through = —^U, can be written as: 


U = -(E sw /ky)exp(x/L E ,v)sin(-k y y/^ lf + 


( 11 ) 


while the related potential energy, W p , reads: 


—^ = a sw — exQ(x / L E v )sin(-k y y + 
mc z c 


( 12 ) 


A plot of the electrostatic potential as a function of space in the surface wave reference 


frame is shown in fig. 10 for the case Rl — 4.3 ( 7 ^ = 4.54, uj/uj pe = 0.22): clearly the initial 
conditions and in particular the initial phase seen by the electron in this potential strongly 
affect the possibility of gaining energy. 
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FIG. 10: Electrostatic potential [Eq. ( 11 )] as a function of space for the case Rl = 4.3 ( 7 ^ = 4.54, 
uj/ uj pe = 0 . 22 , (j) — 0 ). 



21 








If we consider an electron at rest in the laboratory frame, in the boosted frame it will 
have a velocity v y = —v^ and a normalized kinetic energy Wk/ {rnx :2 ) = 7—1 = — 1 . 

If the potential energy is comparable to the kinetic energy, a sw v v /c ~ 7 ^, then for an 
appropriate choice of its initial phase, the electron can be efficiently slowed down in the y 
direction in the boosted frame (and thus accelerated in the laboratory frame). Meanwhile 
it can also gain some perpendicular momentum by exploring the gradient of the potential 
in the perpendicular direction. The net final kinetic energy of the electron after leaving the 
potential has to be of the same order of magnitude of the potential itself: 


W KJ 




m 2 c 2 


+ 


Pf,v 

m 2 c 2 


1 ~ X^iwV^/c, 


(13) 


where % < 2, while we note the final electron momentum components as pf. x and pf )V . It 
should be noted that the limiting value % = 2 corresponds to the case where the electron has 
fully explored the potential in the y direction. It can never be reached in a two-dimensional 
situation due the x component of the SPW field. 

In the laboratory frame, the total energy of the electron is given by: 


with: 


Wtot 
me 2 


7 V 


/ c T 1 T 


v vPf,y \ 
me 2 ) 


(14) 


Px Pxi (15) 

Py = ( Py + mv ip (xa sw v ip /c+ 1)). (16) 


Hence in the limit a sw > 1 and v v ~ c we have: 


Wtot 

me 2 


7 V X 



(17) 


The final energy of the electron in the laboratory frame thus depends on the value and sign of 
pf tV . The best situation would be obtained if Pf, y > 0 and as large as possible, such that the 
electron fully explores the potential in the y direction, provided the potential barrier is high 
enough. In particular, in the limiting case where all the energy is in the parallel direction, 
we will have ^7 ~ a sw and the final energy will be equal to the well-known result for the 
energy gain in a sinusoidal wave [29J 4y v a su , (note that here we have a sw = eE sw /mcaj, where 
E sw is the x component of the field and E y v ~ p v E :iyv ). However, as we already mentioned, 
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TABLE II: Parameters for the different simulations discussed in the paper, pf/mc is the nor¬ 
malised final total momentum obtained for the most energetic electron in the simulation (see figjr]) . 


Uj/0Jpe 

&SW 

r )(fO J sw 

Pf/m e c 

0.22 

8.6 

40 

52 


20. 

91 

130 

0.05 

8.6 

172 

91 


20. 

400 

250 


this limiting situation cannot be reached due to the two-dimensional nature of the potential, 
which will make the electron move away from the surface in the perpendicular direction. It 
is worth noting that the larger the phase velocity, the longer the time the electron can spend 
accelerating within the parallel field before overcoming the wave. Consequently, the effect 
of the finite lifetime of the wave will be important in the higher phase velocity case. An 
upper limit for the time it takes to overcome the wave can be estimated as[28] t,ioj = 2 ir^. 
Thus t sw should be comparable to 77 . 

We can rewrite Eq.(Jl7|) as: 

1 / 1 / ms^rn 

(18) 


Wtot / 

' X X^p^swi 


me * 


where depending on the value and sign of / can be larger or smaller than x- As 
shown in fig| 7 }>c, the perpendicular momentum for the electrons with the highest energies 
is roughly given by p x /mc ~ a sw . Thus going back to the boosted frame and considering 
energy conservation we obtain x ~ 1. If p y > 0, we have x' > 1 which corresponds to the 
case oj/ujp e = 0.22. In this case the numerical value for the normalized final total momentum 
of the most energetic electrons is higher than 7 v a sw , as seen in table II. If instead p y < 0, 
we have %' < 1 which corresponds to the case oj / cu pe = 0.05. In this case, the numerical 
value for the normalized final total momentum of the most energetic electrons is smaller 
than 7 v a sw , as can be seen in table II. 

In order to verify the scaling of the energy gain with a sw we now consider the simulation 
for a sw = 20. The final electron momentum for this case and the three situations previously 
studied is reported in figJTl} As seen in table II, the scaling is confirmed: the normalized 
final total momentum for the most energetic electrons at oj / u pe = 0.22 gives a value for x' = 
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1.4 (instead of 1.3 in the case a sw = 8 . 6 ). As previously pointed out, the finite lifetime of the 
surface plasma wave also affects the value of the normalized final total electron momentum, 
as TdUJ is much larger than t sw uj for u / oj pe = 0.05. As a matter of fact, increasing the lifetime 
by a factor 2 for this case enhances the normalized final total momentum up to pf/mc = 313 
(X' =0.78). 

Finally, a comparison for the electrostatic case is also given in fig(lT]and as expected the 
final momentum is here much less than in the other cases. 
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FIG. 11: Comparison of the final electron momentum in the case a sw = 20 from the ID and 2D 
simulations as a function of the entry phase </> in the surface plasma wave field a) for Rl = 0.9 
( 7 ^ = 1.33, uj/uope = 0.6) b) for Rl — 4.3 ( 7 ^ = 4.54, oo/uj pe — 0.22) and c) for Rl — 20 ( 7 ^ = 20., 
(jo/uo P e = 0.05). o ID simulation, ♦ x compoment of the momentum and * total momentum in the 
2D simulation. Only the electrons in the vacuum side have been considered. 
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CONCLUSION 


In this paper we present different regimes for electron acceleration within the evanescent 
field of a surface plasma wave, which were obtained by means of simulations of the motion 
of test particles. 

For low surface plasma held intensities (characterized by a sw < 1), we show that in the 
electromagnetic regime (characterized by u <C 0 J pe ) the electron motion is essentially ID. 
The electron travels into the vacuum in the direction perpendicular to the surface. This is 
due to the fact that the vacuum evanescence length L E , V of the SPW field is large compared 
to the amplitude of the high-frequency motion of the electron in the wave, v osc /ui. If the 
SPW lifetime is larger than L E}V /i) osc , the electrons with a favorable initial phase will gain 
more energy than provided by the ponderomotive contribution. In the electrostatic regime 
where ao ~ u> pe - \/2, the E y SPW field components induce positive 2 D effects, which amplify 
the electron energy gain. 

In the regime where relativistic effects are no longer negligible (for a sw >1), the relative 
weight of the perpendicular versus parallel SPW field components may limit the interaction 
of the electron with the full two-dimensional field of the surface wave. Nevertheless, we find 
optimal conditions for acceleration in which the electron can still gain a large amount of 
energy, mainly by acceleration parallel to the surface. This is achieved when the electrons 
acquire a parallel velocity close to the phase velocity, v^, of the wave propagating along 
the surface. This is observed in the so-called relativistic electromagnetic regime where 
7^ = V \/l — (w/ c ) 2 is large and uo -C co pe . Reaching this optimal regime also depends on 
the lifetime of the surface plasma wave. 

In this regime an electron initially at rest is self-injected and phase-locks on the vacuum 
side. A significant contribution to the self-injection is provided by the y component of the 
vxB force (v x B) resulting in favorable 2 D effects in contrast to a purely electrostatic wave. 
This way the electron gains a substantial amount of parallel energy. In the boosted frame 
on the vacuum side the B field becomes equal to zero while E y v is invariant. However, even 
if in this frame the particle is not initially in phase with the wave and j lfi > a sw v v /c, the 
electron can still gain a large amount of energy. This happens while it is subjected to the 
parallel field during the accelerating part where it is rolling away from the surface before 
the field changes its sign and becomes decelerating. This confirms the favorable effect of 
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two-dimensional dynamics. On the plasma side (this case has not been presented here) this 
mechanism is less efficient because the perpendicular component of the SPW electric field is 
smaller than in vacuum. Also its spatial extension is smaller. Moreover, on the plasma side 
the problem does not become electrostatic in the boosted frame. 

We point out that an optimum energy gain could also be reached for external injected 
electrons with an high initial velocity and an appropriate angle of injection. 

Finally, we want to outline that we have examined only electrons that are accelerated and 
are spending their dynamical history on the vacuum side. Depending on the phase, some 
electrons can also spend their dynamical history inside the plasma. Since the electric field 
perpendicular to the surface is discontinuous and changes sign we found that some electrons 
can oscillate back and forth from the surface to the plasma. After some oscillations they 
may stay on one side or another. However, we decided not to study such motion in the 
present work, since the numerical study cannot be precise in the presence of a discontinuity. 
Moreover, a more realistic situation including plasma temperature effects will easily smooth 
such a discontinuity, without otherwise significantly changing the surface wave solutions [3D] • 
Other important effects that are missing in this simple model, but have been observed - e.g. 
in PIC simulations [30] - are the presence of a charge space field at the plasma surface left 
by the electrons ejected from the surface, or the presence of a quasi-static magnetic field 
PSl El]. Considering these effects is beyond the scope of this paper and will be the subject 
of further studies. 

One of ns (C.R.) would like to acknowledge fruitful discussions with F. Amiranoff, A. 
Macchi, P. Mora and D. Pesme. 
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